Quantum quenches in disordered systems: 
Approach to thermal equilibrium without a typical relaxation time 
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We study spectral properties and the dynamics after a quench of one-dimensional spinless fermions with 
short-range interactions and long-range random hopping. We show that a sufficiently fast decay of the hopping 
term promotes localization effects at finite temperature, which prevents thermalization even if the classical 
motion is chaotic. For slower decays, we find that thermalization does occur. However, within this model, the 
latter regime falls in an unexpected universality class, namely, observables exhibit a power-law (as opposed to 
an exponential) approach to their thermal expectation values. 
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The study of the dynamics and conditions for thermaliza- 
tion in isolated many-body systems has a long history in the- 
oretical physics 01]. In classical physics, the requirements for 
thermalization are well understood. Boltzmann's ergodic hy- 
pothesis holds only if the motion of the individual particles is 
fully chaotic. The situation in quantum systems is less clear. 
Quantum time evolution is described by a linear differential 
equation so that dynamical chaos does not occur. The de- 
velopment of the theory of quantum chaos in the 1980's and 
1990's brought a new language and techniques to tackle this 
problem 10]. With respect to thermalization, this effort crys- 
tallized in two main results: the so-called Berry's hypothesis 
Jit], which states that eigenstates of a classically chaotic sys- 
tem can be written as a finite sum of plane waves with ran- 
dom coefficients, and the eigenstate thermalization hypothe- 
sis (ETH) proposed by Deutsch 0] and Srednicki fU, which 
states that observables in individual eigenstates of a generic 
many-body system already exhibit thermal behavior. 

Until a few years ago, technical difficulties prevented a sys- 
tematic study of the proposals above. However, recent ad- 
vances in cold gases systems, together with the enhancement 
of computer power and the development of novel computa- 
tional methods, are making possible a more quantitative com- 
parison. This has dramatically boosted the theoretical and 
experimental interest in non-equilibrium dynamics in general 
and thermalization in particular @. For instance, in Ref. 01, 
it was shown experimentally that, after a quench, the momen- 
tum distribution of a gas of bosons trapped in a quasi-one- 
dimensional (ID) geometry did not relax to the thermal pre- 
diction. Although this is expected in an integrable system JUl, 
it was surprising that such an effect could be seen experimen- 
tally. The ETH, on the other hand, has been confirmed nu- 
merically for nonintegrable systems |9J, and has been shown 
to break down when approaching integrable points HoHUl. 

A better understanding of the conditions for thermalization 
would not only put the foundation of quantum statistical me- 
chanics on a more solid ground but also have a strong im- 
pact in different fields, from cold gases to cosmology, where 
non-equilibrium dynamics play a key role. Here, we provide 



further insights on this problem. We show that many-body 
localization effects foil invalidate the expectation that classi- 
cal chaos leads to thermalization of the quantum counterpart. 
(For recent connections between the effect of localization in 
Fock space and thermalization in spin systems, see Ref. fTill .) 
We also put forward a route toward thermalization in quan- 
tum systems characterized by power-law relaxation dynamics. 
We support these results by numerical calculations in the fol- 
lowing ID spinless fermions system, with long-range hopping 
and short-range interactions, 
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where fj creates a fermion in the site j, and h } is the number 
operator in the site j. The hopping term is built by means of a 
Gaussian random distribution, with zero mean = 0, and 
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In the limit V = 0, the properties of (Q3 depend on a but 
not on /3 > lfl5ll . For a < 1, eigenstates are delocalized 
and the spectral correlations are described by Wigner-Dyson 
(WD) statistics. For a > 1, eigenstates are localized and spec- 
tral correlations are described by Poisson statistics. For a = 1, 
eigenstates are multifractal and spectral correlations are in- 
termediate between WD and Poisson iflril [nil . We then fix 
j5 = 0. 1 and V = 1. (The latter sets the unit of energy through- 
out this Rapid Communication, and £g = 1.) These values 
were chosen to minimize finite-size effects (j3 <C 1) and at 
the same time to avoid the trivial insulator limit that occurs if 
the potential energy is much larger than the kinetic one ifTTll . 
Likewise, if the interaction is too weak, thermalization may 
not occur ifioll . Finally, we note that these types of long-range 
Hamiltonians have been used to model systems with strong 
dipolar interactions lfl8ll . 

We first use a finite-size scaling analysis to investigate the 
effect of the many -body interactions on the spectrum lfl7ll . 
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FIG. 1. (Color online) Scaling variable 77 [see Eq. (f3]> ] as a function 
of a for different system sizes but the same filling factor, 1 /3. For 
a ~ 1 (°t 1-2), ?7 decreases (increases) with the system size. This 
is a signature of a metallic (insulating) phase. The number of realiza- 
tions is 10000, 10000, and 400 for L = 12, 15, and 18, respectively. 

This is a powerful tool to study many-body localization in the 
presence of interactions ifTTl [l3ll . We compute the eigenval- 
ues of the Hamiltonian ([1) for different sizes and values of 
a utilizing standard diagonalization techniques. In all cases, 
the filling factor p/L = 1 /3, where L is the system size and 
p is the number of particles. The spectrum thus obtained is 
appropriately unfolded, i.e., it is rescaled so that the spectral 
density on a spectral window comprising several level spac- 
ings is unity. The number of disorder realizations considered 
is such that statistical fluctuations are negligible. As a scaling 
variable we choose the function 77 11711 related to the variance 
var = (s 2 ) — (s) 2 of the level spacing distribution P(s). P(s) 
is the probability of finding two neighboring eigenvalues at 
a distance s = (e, + i — £,)/A, where A is the local mean level 
spacing, and 

77 = [var - var WD ] / [ var P - var W D] • (3) 

varwD = 0.286 (varp = 1) is the value of the variance for a 
disordered metal (insulator) in the thermodynamic limit and 
(s n ) = J s"P(s) IU9I1 . An increase (decrease) of 77 with L sig- 
nals that the system will be an insulator (metal) in the thermo- 
dynamic limit. 

Figure Q] depicts results for 77 vs a for different sizes. It is 
apparent that for a > 1 .2 the parameter 77 increases with sys- 
tem size, as is expected of an insulator. Hence, localization 
takes place in the interacting system, in contrast to the classi- 
cal counterpart, which is chaotic for any a. For a < 1, on the 
other hand, 77 behaves as expected of a metal. To be certain 
whether the system is metallic for a « 1 , much larger systems, 
currently not accessible numerically, need to be studied. 

We now investigate the thermalization properties of the 
Hamiltonian ([1]). We aim to (i) identify a region of a's for 
which the system does not thermalize even though the clas- 
sical counterpart does, (ii) see how that region relates to the 
localization regime found by the spectral analysis, (iii) inves- 
tigate the microscopic origin of the lack of thermalization, and 
(iv) study the approach to equilibrium in the region of a's for 
which thermalization eventually occurs. 



In order to proceed, we first note that time scales and finite- 
size effects relevant to the study of quantum thermalization 
may depend on the observable and particle statistics [QUO. 
However, for few-body observables, it is generically expected 
that thermalization occurs away from integrability. Here, we 
report results for two of those observables: the momentum 
distribution function [n(k)] and the density-density structure 
factor [N(k)], 

l,m Ijn 

which are the Fourier transforms of the one-particle and 
density-density correlation matrices ll20ll . respectively. Both 
observables can be measured in ultracold gases experiments. 

In all the cases shown below, we start from an eigen- 
state of the Hamiltonian ([TJ in a certain realization of the 
disorder. Then, we change to another disorder realization 
for the same a, and study the time evolution of the ini- 
tial state with this final Hamiltonian, H^. This procedure 
is usually known as a quench. The initial state [^(O))] is 
selected such that the time evolving system has an energy 
E = ( x P(0)|//fi n | l P(0)}, which, for every disorder realization, 
is the same as the one of a thermal state with temperature 
T = 10.0 (E = Tr{e^fl-/ r i/ fin }/Tr{e^fl-/ :r }). This yields 
energies that fall close to the center of the spectrum of the fi- 
nal Hamiltonian. In what follows, Ojj are the matrix elements 
of a given observable in the eigenstates of the final Hamil- 
tonian, Oij = and C; is the overlap between the 
initial state and Cj = (y/yl^O)). 

In order to determine whether thermalization occurs fol- 
lowing the quench, one needs to find a meaningful quan- 
tity to compare with the microcanonical (thermal) average, 
Omicro = -&^HiOji, where JV^e is the number of states in 
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the microcanonical window (centered around E, and with AE 
selected such that the average is robust against small changes 
of AE). If relaxation takes place for the observables of in- 
terest, and the spectrum is nondegenerate, the infinite time 
average (also known as the diagonal ensemble prediction) 
^diag = L; \Cj\ 2 Ojj is the right choice JH]. We first compute 
the normalized difference between these two ensembles, 

. La: I ^diag (&) ^micro 

LtOdiagO) 

and then average it over different disorder realizations to ob- 
tain (AO)dis- Note that here, and in what follows, by "O" we 
mean "n" or "N" for the momentum distribution and structure 
factor, respectively. 

In Fig. 121 we depict (A«)di s and (AiV)dis for different val- 
ues of a vs system size. Thermalization occurs if (AO)dis 
vanishes in the thermodynamic limit. A nonzero value of 
(A(9)di s in this limit indicates that the observable O relaxes 
to a nonthermal expectation value. For a > 1, a weak size 
dependence is observed for the largest system sizes we can 
study, with (A(9)di s likely saturating to non-zero values for 
a > 1.2. Therefore, thermalization is not expected to occur in 
this regime. Interestingly, as the value of a decreases below 
a ~ 1, the normalized differences exhibit a fast decrease for 
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FIG. 2. (Color online) {A/7)d; s and (AN) fa [see Eq. as a func- 
tion of the system size, for a = 0.6, 0.8, 1.0, 1.2, and 1.4. Points 
correspond to (L,p) = (9,3), (12,4), (15,5), and (18,6). The disor- 
der average is performed over at least 8 , 500 different realizations for 
L < 15, and 1020 for L = 18. The classical dynamics is chaotic but 
there is no thermalization for large a due to localization effects. 



the smallest systems shown. They become much smaller than 
those for a > 1 .2 for the largest system sizes accessible here, 
for which (A<9)di s is very close to zero within our error bars 
and still decreasing with increasing system size. These results 
suggest that thermalization occurs in this regime. 

In order to further support the conclusions of the finite- 
size scaling analysis we look at the actual diagonal and mi- 
crocanonical expectation values of observables for several 
quenches. Results for n(k = 0) are shown in Fig.|3]as a func- 
tion of the energy. In all regimes, the microcanonical results 
can be seen to be almost independent of the energy, while the 
diagonal ensemble results exhibit fluctuations that increase as 
a increases. Hence, increasing a increases the difference be- 
tween the infinite time average and the microcanonical results, 
as well as increases the dependence of the infinite time aver- 
age on the initial state selected. 

A natural question that follows is whether the absence of 
thermalization, as well as the sensitivity to the initial state se- 
lected, for large a, is related to the breakdown of ETH (as seen 
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FIG. 3. (Color online) Microcanonical and diagonal results for 
n(k = 0) in 51 different quenches. The systems have 18 sites and 
six particles, with a = 0.6, 1.0, and 1.4. 
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FIG. 4. (Color online) (A/i^ ax ) dis and (AA^ ax ) dis (main panels), and 
(Aii*J°}di s and {AA^ vs )d; s (insets) vs system size. Lines are the same 
as in Fig. [2] As a consequence of localization effects, ETH does not 
hold for large values of a. 

in clean systems approaching an integrable point iflQjfTTll ) or 
it is rather related to some atypical properties of the overlaps 
Cj. To answer that question, we compute the normalized dif- 
ference between the observable in each eigenstate and in the 
microcanonical ensemble, 

AOii = ywi 7T\ ' W 

This allows us to determine, for each disorder realiza- 
tion, the maximal difference within the microcanonical win- 
dow A(9™ x = Max[AO,-,]A£ as well as the average AO™ s = 
-p— y\ AOu ■. We then average both quantities over different 
disorder realizations to obtain (A<9™ ax )d; s and (AO^ vs )dis- 

In the main panels in Fig. @] we depict (An™ ax )di s and 
(AJV™ ax )dis vs L for different values of a. ETH holds when 
(A(9™ ax )di s — > for L — >• oo. In the range of sizes that we can 
study, this behavior is apparent for a < 1 . For a > 1 .2, we find 
clear indications that ETH fails, which can be understood as 
a result of localization induced by disorder |[T3l . A very sim- 
ilar behavior is observed in the insets of Fig. |4] which show 
(A«^ yg )di s and (AA^ vs )di s . In the region a w 1, on the other 
hand, our results are not conclusive. This is an interesting 
problem for future works as, in the noninteracting limit, a = 1 
corresponds to a metal-insulator transition characterized by 
multifractal eigenstates. We speculate that fluctuations at all 
scales associated with multifractality may lead to interesting 
behavior in the many-body properties of the system. 

The robustness of the results for (A(9Jf ax )di s and (A<9^ vg )di s , 
as well as their clear correlation with the thermalization indi- 
cators in Fig. |2 allow us to conclude that: (i) the lack (oc- 
currence) of thermalization is directly related to the failure 
(validity) of ETH, and (ii) that ETH holds and thermalization 
occurs only for values of a < 1 . For a greater than, and not 
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FIG. 5. (Color online) Time evolution of (8n{t))^ s and (SN(t))^ s 
[see Eq. 0] for different a's. Thick lines are for 18 sites and six 
particles. Thin solid lines are power-law fits to the data. Other thin 
lines are the corresponding results for 15 sites and five particles. An 
average over 8500 (1020) realizations has been carried out for the 
15-site (18-site) system. 



too close to, one, the quantum system will not thermalize even 
though the dynamic of the classical counterpart is chaotic. 

A fundamental problem that has not been addressed in pre- 
vious computational studies-due to large fluctuations that oc- 
cur during the time evolution in exact diagonalization stud- 
ies (because of finite-size effects) and to the limited times 
accessible in time-dependent density-matrix renormalization 
group (t-DMRG) approaches Il2lll -is that of how observables 
approach their thermal values during the relaxation dynam- 
ics. The naive expectation is that the approach should be ex- 
ponential as in classical systems, where a few collisions per 
particle suffice for the system to relax to thermal equilibrium. 
However, to the best of our knowledge, such a behavior has 
yet to be seen in the experiments or numerical simulations of 
isolated systems in the quantum regime. Disordered systems 
provide a unique opportunity to address this problem as the 
average over disorder realizations reduces dramatically fl uctu- 
ations due to finite-size effects. In what follows, we compute 
the time evolution of the difference 

ao( , )= M-y>i m 



and then average it over different disorder realizations to ob- 
tain (50(f))dis OH] . In Fig. [5] we show results for (5«(f))dis 
and (5iV(f ))di s as a function of time, t, for three different val- 
ues of a and the two largest system sizes that we can study. 
For a > 1, an extremely slow relaxation dynamics can be 
observed, and the system may never reach the infinite time- 
average prediction in any experimentally relevant time scale. 
For a < 1, the relaxation dynamics seen in those plots is 
quite unexpected. We find that (<5n(f))di s and (&/V(f ))di s ex- 
hibit a clear power-law decay (°= t ~f) The region over which 
the power-law decay is observed extends over a decade and 
increases with increasing system size. As the value of a 
decreases (and localization effects decrease) the exponent y 
of the power law increases. However, no typical time-scale 
emerges during relaxation ll22ll . This indicates an unexpected 
route to thermal equilibrium in many-body quantum systems 
characterized by a power law rather than an exponential decay. 

After these results, it is natural to speculate whether such 
power-law behavior also occurs in clean systems. Theoreti- 
cally, it is well known that in the semiclassical limit classi- 
cal cantori l23ll . remnants of the Kolmogorov-Arnold-Moser 
(KAM) tori induce slow diffusion in phase space and power- 
law localization of eigenstates in the one-body problem l24ll . 
Therefore, it is plausible that in certain region of parameters 
the approach to equilibrium in systems controlled by cantori 
is also power-law-like. Interestingly, indications of power-law 
relaxation have already been seen in classical systems rf25ll 
and, experimentally, in a strongly correlated one-dimensional 
Bose gas I2r3l . 

In conclusion, we have studied an interacting many-body 
disordered system that exhibits a transition between metallic 
and insulating behavior. Remarkably, we have identified a re- 
gion of parameters (a > 1.2) in which, due to localization ef- 
fects, ETH fails and thermalization does not take place even if 
the system is nonintegrable 111 411 . For a < 1, ETH is valid and 
thermalization occurs. Furthermore, in the latter regime, we 
have found a route toward thermal equilibrium characterized 
by a power-law approach to the thermal expectation values 
and, hence, by the lack of a well-defined equilibration time. 
The relevance of this scenario to experiments with ultracold 
gases, as well as clean strongly correlated systems, is a topic 
that requires future exploration. 
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